Error-resistant Single Qubit Gates with Trapped Ions 
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Coherent operations constitutive for the implementation of single and multi-qubit quantum gates 
with trapped ions are demonstrated that are robust against variations in experimental parameters 
and intrinsically indeterministic system parameters. In particular, pulses developed using optimal 
control theory are demonstrated for the first time with trapped ions. Their performance as a function 
of error parameters is systematically investigated and compared to composite pulses. 
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In order to experimentally implement a device capable 
of performing fault-tolerant universal quantum compu- 
tation (QC), quantum gate operations involving one or 
multiple qubits have to be carried out with demandingly 
high accuracy (see, for instance, [EH])- According to 
recent theoretical investigations, the experimentally re- 
quired accuracy of quantum gates for fault-tolerant uni- 
versal quantum computation no longer seems daunting 
or even prohibitive [2f- But still, the desired error prob- 
ability per gate (EPG) should be as small as possible 
in order to keep the experimental overhead necessary for 
quantum computation within a feasible limit. Thus a low 
error probability is prerequisite for scalable fault-tolerant 
QC. 

Any quantum algorithm can be decomposed into a se- 
quence of unitary operations applied to individual qubits 
(single-qubit gate) and conditional quantum dynamics 
with at least two qubits Multi-qubit gates (involv- 
ing two or more qubits) are synthesized by applying a 
sequence of elementary unitary operations on a collec- 
tion of qubits. Each of these elementary operations is 
often similar, or identical, to what is needed for single- 
qubit gates, and therefore each operation has to be imple- 
mented with an error probability well below the tolerable 
EPG characterizing the full gate operation. 

If electrodynamically trapped ions are used as qubits, 
then a unitary operation amounts to letting ions inter- 
act with electromagnetic radiation with prescribed fre- 
quency, phase, amplitude, and duration of interaction in 
order to implement quantum gates. Recently, impressive 
experimental progress was demonstrated in entangling 
up to eight ions, and performing 2-qubit quantum gates 
0, H, El- Architectures allowing for scalable QC with 
trapped ions have been proposed (e.g., [7[), and build- 
ing blocks necessary for achieving this ambitious goal are 
currently being investigated using various types of ions. 

The error budget, for instance, of the geometrical 
phase gate demonstrated in Q is dominated by the fre- 
quency and amplitude uncertainty of the laser light field. 
These errors are also responsible for a part of the EPG 
of the controlled-NOT gate reported in Q. If an "ion 
spin molecule", that is, trapped ions coupled via a long 
range spin-spin interaction, is to be used for quantum in- 



formation processing, then the exact transition frequency 
of a particular ionic qubit depends on the internal state 
of other ions [8|. Therefore, here too, it is important to 
have quantum gates at hand that are insensitive to the 
detuning of the radiation driving the qubit transition. 

Here, we demonstrate single qubit gates with trapped 
ions that are robust against experimental imperfections 
over a wide range of parameters. In particular it is shown 
that errors caused by an inaccurate setting of either fre- 
quency, amplitude, or duration of the driving field, or 
of a combination of these errors are tolerable (in terms 
of a desired accuracy of quantum gates) when a suit- 
able sequence of radiation pulses is applied instead of, 
for instance, a single rectangular 7r-pulse. Thus an essen- 
tial prerequisite for scalable quantum computation with 
trapped ions is demonstrated. 

We shall show results of using shaped pulses which 
were developed using optimal control theory. Optimal 
control theory relies on a generalization of the classical 
Euler-Lagrange formalism. Using discrete rather than 
continuous trajectories in the configuration space of the 
quantum mechanical system, and constraining the tra- 
jectories to satisfy the Bloch equations, the cost function 
to be optimized, is the transfer efficiency over a range of 
frequency detunings between qubit and applied radiation 
and amplitudes of the radiation 0] . In addition, compos- 
ite pulses are realized here that are specifically designed 
to tackle off-resonance errors, or designed to tackle pulse 
length errors or power discrepancies of the driving field 

m 

The two level quantum mechanical system used as 
qubit is realized on the |0) = Si/ 2 (F=0) <-» Si/a(F=l, 
mj?=0) = |1) transition in a single 171 Yb + ion with 
Bohr frequency u>o confined in a miniature Paul trap 
(diameter of 2 mm), driven by microwave radiation 
with frequency to close to 27r x 12.6 GHz and Rabi fre- 
quency w 2tt x 10 kHz. The time evolution of the 
qubit is virtually free of decoherence, that is, transver- 
sal and longitudinal relaxation rates are negligible, and 
is determined in a rotating frame after the rotating 
wave approximation, by the semi-classical Hamiltonian 
H = \bo z + |fi (cr+e- 4 * + cr-e^) . Here, cr± are the 
atomic raising/lowering operators, a z is a Pauli matrix, 
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and 5 = luq — u> is the detuning of the applied radiation 
with respect to the atomic transition. Imperfect prepa- 
ration of the |0) state by optical pumping, limited here 
the purity of the initially prepared state, such that the 
initial density matrix (before coherent interaction with 
microwave radiation) is given by pi = <xi|l)(l| + ao|0)(0| 
with typical values a\ = 0.1 and ao = 0.9. 

The ion is produced from its neutral precursor by 
photoionization using a diode laser operating near 399 
ran. Laser light near 369 nm driving resonantly the 
S1/2 F=l <-> P1/2 F=0 transition in Yb + is supplied by 
a frequency doubled Ti:Sa laser, and serves for cool- 
ing and state selectively detecting the ion. A diode 
laser delivers light near 935 nm and drives the D 3 / 2 - 
[3/2] ^transition to avoid optical pumping into the 
metastable D 3 / 2 state du ring the cooling and detection 



periods (for details see 12|). For quantifying errors 



in detuning we will use the scaled detuning, / = 5/Q 
whereas errors in pulse area shall be represented by 
g = A8/9, where A9 = 9' — 9 with the desired pulse 
area, 9 = L £ldt and the actual pulse area, 9' when T 
or fl are not set perfectly. The fidelity of the qubit state 
\6 m , <pm) = cos (0 m /2) | 0) + e^™) sin (9 m /2) | 1) that is 
obtained after applying a microwave pulse is given by 
F =| {9,4> | 9 m ,4> m ) | 2 with \9,(p) being the state that 
would be obtained, if the pulse were perfect. Thus, for a 
9 = 7r-pulse the fidelity is given by F„ =\ sin(# m /2) | 2 . 
Impure initial preparation (ai > 0) limits the maxi- 
mum fidelity that can be obtained with a 7r-pulse to 
F-n = ao — a\ . In order to determine the fidelity of the 
state obtained after a gate that is supposed to leave the 
qubit with a well defined phase <fi {e.g. a 7r/2-pulse), 
the phase <f> m needs to be determined in addition to 
9 m . A Ramsey- type experiment allows for measuring 
both angles: First an ideal Ramsey sequence is car- 
ried out (i.e. two successive ideal 7r/2-pulses with vary- 
ing phase <E> of the second pulse) yielding interference 
fringes in the population of state |1), a\ as a function 
of $. Then, the first 7r/2-pulse is replaced by a possi- 
bly non-ideal pulse sequence leaving the qubit in state 
\9 m ,4>m), and again interference fringes, are recorded. 
Now the population ai($) detected in state |1) is given 
by ai(<f>) = 1/2(1 + sin(# m ) cos(</> m - $)), and from a fit 
of the data one obtains 9 m and <p m . 

The basic measurement sequence (labeled sequence A) 
for determining the fidelity as a function of / and g of 
a shaped pulse that ideally gives a rotation with 9 = n 
is as follows: i) A single ion is prepared in the |0) state 
by optical pumping through illumination with 369 nm 
light for 20ms. ii) A microwave shaped pulse with con- 
trolled error, that is, known values of / and g, is applied, 
iii) Again, the ion is illuminated for 4 ms with 369 nm 
light for state selective detection, iv) The ion is laser 
cooled by applying microwave and laser radiation simul- 
taneously. This sequence comprising steps i) through iv) 
is then repeated (sequence B), except that in step ii), 
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FIG. 1: Upper graphs: Ramsey interference fringes from se- 
quence A' (left hand side) and sequence B' (rhs, / = — 1, 
g = 0). A fit (indicated by a solid line) of the experimental 
data (crosses) gives m = 3.92 ± 0.19, m = 0.418 ± 0.34 (see 
text). Lower: Relative amplitude (lhs) and phase (rhs) of a 
pulse consisting of 645 individual steps evolving over time. 
This pulse was developed using optimal control theory [j~H . 



for direct comparison, the shaped pulse is replaced by a 
rectangular pulse that would give 9 = 7r, if / = = g. 
Then, i) through iv) is repeated again (sequence C) with 
an ideal (i.e., f = = g) n/2 rectangular pulse to control 
state selective detection in each measurement cycle. Fi- 
nally, i) through iv) is repeated a fourth time (sequence 
D) leaving out step ii) in order to monitor the initial 
preparation in terms of the coefficient a Q . Typically, the 
full procedure (sequences A through D) outlined here is 
repeated 700 times for a given pair of / and g values. 

When measuring the performance of shaped pulses and 
pulse sequences that ideally yield a rotation of 8 = n/2 
and, for instance, <f> = — n/2, then the basic sequence, A' 
is as described above, except that in step ii) two ideal 
(/ = = g) n/2 pulses are applied, the second one with 
variable phase $ (this yields Ramsey fringes as a refer- 
ence). Then, sequence B' is performed by replacing the 
first 7r/2 pulse in step ii) by a shaped pulse with con- 
trolled error. Sequence C is obtained by replacing the 
first 7r/2 pulse in ii) with a rectangular pulse subject to 
the same errors as the shaped pulse in B'. Sequences A' 
through C are repeated 20 times while increasing the 
value of $ by 27r/20. Then sequences C and D are ap- 
pended and the complete procedure is repeated 50 times. 
As an example, Fig. [1] upper panels illustrates how an- 
gles 9 m , 4> m were determined for a specific pair of / and 
g values. 

A completed measurement returns two grids of fideli- 
ties with the points on each grid defined by different val- 
ues of / and g. One grid corresponds to a simple pulse 
(obtained from sequences B or C), that is, the amplitude 
as a function of time has a rectangular shape with the 
desired pulse area, 9 = fit. This pulse results in a perfect 
rotation by an angle 9 only for / = and g = 0. The 
second grid of fidelities corresponds to a shaped or com- 
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FIG. 2: Upper panel: Experimental fidelity for the shaped 
7r/2-pulse (His) shown in Fig. [1] and of a rectangular n/2- 
pulse (rhs) as a function of detuning / and amplitude g. The 
average statistical error of the measured points for the shaped 
pulse is Gav = 0.067, whereas a av = 0.086 for the rectangular 
pulse. Lower panel: A white rectangle indicates F/F m > 0.90 
for the shaped and the rectangular pulse, respectively. 

posite pulse (obtained from sequence A or B'). In the ex- 
perimental results shown below, crossing points between 
gridlines represent measurement points. The shaded ar- 
eas are obtained by linear interpolation between points. 

The pulses developed using Optimal Control Theory 
and demonstrated here were designed for off-resonant 
errors described by —1 < / < 1, power variations up 
to ±40%, Rabi frequency SI = 2tt x 10 kHz, and are 
made of 0.5/is steps. In each step, in general, ampli- 
tude and phase of the radiation is changed. As an ex- 
ample the time evolution of phase and amplitude of a 
shaped pulse consisting of 645 individual steps is shown 
in Fig. [T] [ll[ . This pulse is supposed to yield a rotation 
of 9 = it/2, (f> = — 7r/2 that is robust against variations in 
/ and g. 

The experimentally determined fidelity of this shaped 
pulse as a function of / and g is displayed in Fig. For 
reference, Fig. [2] also shows the experimental fidelity ob- 
tained from a simple rectangular tt/2 gate for the same 
parameter range. For / = = g the simple pulse yields 
the maximum fidelity F m — 0.899 ± 0.065, and the range 
of / and g for which F/F m > 0.90 is indicated in Fig. 
[2] lower panel by white rectangles. The same maximum 
fidelity, F m = 0.900 ± 0.064, is also reached using the 
shaped pulse. However, this maximum value is main- 
tained over a much wider parameter range as is evident 
from Fig. [2J thus demonstrating the robustness of the 
shaped pulse against experimental errors and intrinsic 
imperfections. 

Fig. [3] displays the experimental fidelity obtained by 
using a shaped 7r-pulse consisting of 445 steps with vari- 
able phase and amplitude subject to controlled errors. 
Again, a rectangular pulse serves as an experimental ref- 
erence giving a maximum fidelity F m = 0.838 ± 0.030 for 
/ = = g that rapidly decreases for increasing |/| or \g\. 
The shaped pulse, in contrast, maintains the maximal 



FIG. 3: Upper panel: Experimental fidelity of a shaped tt- 
pulse (lhs, 445 steps, a av = 0.028), and the corresponding 
rectangular 7r-pulse (rhs) as reference (a a v = 0.020). Lower 
panel: White rectangles show F/F m > 0.96. 



possible fidelity over a wide range of parameters as is ev- 
ident in Fig. [3l A shaped 7r-pulse consisting of 835 steps 
was also implemented resulting in a parameter area indi- 
cating robustness against errors that extends even further 
(not shown). 

If one uses shaped pulses obtained from optimal control 
theory, even with parameters that considerably deviate 
from their ideal values, one still attains accurate perfor- 
mance of single qubit gates. Now, we compare the per- 
formance of these pulses to some composite pulses that 
have been devised to be effective against off-resonance, 
amplitude, and pulse length errors. The resulting fidelity 
grids were measured by sequentially performing the com- 
posite case and the simple case for each point on the grid, 
and then repeating 300 times. In terms of the previously 
used description, sequence A and sequence B are per- 
formed for each point on the grid, they are followed by 
C and D, then the complete measurement procedure is 
performed 300 times. 

The composite pulse of type CORPSE(Compensation 
for Off-Resonance with a Pulse Sequence) was derived 
with the aim of combatting off- resonant errors [rol |. It 
consists of three pulses where the first and the third 
pulses have equal phase values $ and the phase of the sec- 
ond pulse differs by w. The nutation angles of the three 
pulses making up a it composite pulse are 420°, 300°, 60°. 
Indeed, as shown in Fig. [4] upper panel the CORPSE 
7r-pulse extends the experimentally determined range of 
detuning / over which a fidelity F/F m > 0.96 (F m = 
0.915 ± 0.043) is maintained as compared to the rectan- 
gular pulse shown on the right of Fig. [U However, for 
compensating pulse area errors this pulse sequence is less 
effective as is also evident from Fig. @] 

The SCROFOLOUS (Short Composite Rotation For 
Undoing Length Over and Under Shoot) pulse was de- 
rived with the aim of compensating for pulse area errors. 
These pulses were derived by restricting the phase and 
and nutation angles such that $1 = $3;0i = 6*3 with O2 
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FIG. 4: Upper panel: Measured fidelity for a CORPSE-type 
composite pulse (lhs, 8 = n, a av = 0.033) with the corre- 
sponding rectangular pulse for reference (rhs, a a v = 0.030). 
Also shown are measured grids in terms of F/F m > 0.96 where 
white is true and black is false. Lower panel: Measured fi- 
delity of a SCROFOLOUS-type composite pulse (lhs, 8 = n, 
Oav = 0.0327) with the corresponding rectangular pulse for 
reference (rhs, a a v = 0.0299). White squares in the graphs 
below indicate F/F m > 0.96, F m = 0.915 ± 0.043. 



and $2 unbound. The angles are chosen such that once 
again the fidelity is least effected by the presence of er- 
rors [nj. For a 7r-pulse these angles are 6\ = 180°, 9 2 = 
180°; $i = 60°, $ 2 = 300°. We show in Fig. H lower 
panel that errors in pulse area caused by power fluctua- 
tions are well compensated for. On the other hand, de- 
tuning errors have the same detrimental effect as is the 
case with a rectangular pulse. Composite 7r pulses, with g 
describing time errors rather than amplitude errors were 
also implemented (not shown). 

The BBI (Broadband) composite pulse comprises a se- 
quence W = I80$i — 360^2 — 180$3. When the desired 
rotation is R(#), the complete rotation can be performed 
as R{9) -W,W- R{9) or (9/2) - W - R(6/2). The mea- 
surement procedure follows the procedure outlined above 
for the shaped 7r/2-pulse. Shown in Fig. [5] are the results 
of using the BBI composite pulse R(9/2) - W - R{6/2) 
with 9 = 7r/2. The experimental data show that BBI is 
effective for compensation of errors in both detuning and 
pulse area. The parameter range over which error resis- 
tant pulses are effective is greater than that of a simple 
pulse but is more restricted than with the shaped pulse 




FIG. 5: Upper panel: Experimental Fidelity for the BBI (§) — 
W - i?(f ) sequence (lhs), where 8 = f (left, a av = 0.072) 
and for a corresponding rectangular pulse (rhs) . Lower panel: 
Measured grids in terms of F/F m > 0.90, F m = 0.936 ± 0.052 



shown in Fig. [2] 

A comparison of the performance of shaped pulses de- 
veloped using optimal control theory with simple rectan- 
gular pulses or composite pulses reveals an evident ad- 
vantage of these shaped pulses in terms of robustness 
against experimental errors and indeterministic system 
parameters, while the lengths of both types of pulses are 
comparable (here, with il = 2ir x 10 kHz, for instance 
the pulse in Fig. [3] takes 223 /is compared to 217 /xs for a 
CORPSE pulse of Fig. H]). This will make shaped pulses 
based on optimal control theory an important tool in or- 
der to achieve quantum gates with trapped ions with low 
error probability and thus come a step closer to fault- 
tolerant quantum computing. 
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